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ABSTRACT 



Context. Previous molecular gas observations at arcsecond-scale resolution of the Seyfert 2 galaxy M51 suggest the 
presence of a dense circumnuclear rotating disk, which may be the reservoir for fueling the active nucleus and obscures 
it from direct view in the optical. However, our recent interferometric CO(3-2) observations show a hint of a velocity 
gradient perpendicular to the rotating disk, which suggests a more complex structure than previously thought. 
Aims. To image the putative circumnuclear molecular gas disk at sub-arcsecond resolution to better understand both 
the spatial distribution and kinematics of the molecular gas. 

Methods. We carried out CO(2-l) and CO(1-0) line observations of the nuclear region of M51 with the new A configu- 
ration of the IRAM Plateau de Bure Interferometer, yielding a spatial resolution lower than 15 pc. 
Results. The high resolution images show no clear evidence of a disk, aligned nearly east-west and perpendicular to the 
radio jet axis, as suggested by previous observations, but show two separate features located on the eastern and western 
sides of the nucleus. The western feature shows an elongated structure along the jet and a good velocity correspondence 
with optical emission lines associated with the jet, suggesting that this feature is a jet-entrained gas. The eastern 
feature is elongated nearly east-west ending around the nucleus. A velocity gradient appears in the same direction with 
increasingly blueshifted velocities near the nucleus. This velocity gradient is in the opposite sense of that previously 
inferred for the putative circumnuclear disk. Possible explanations for the observed molecular gas distribution and 
kinematics are that a rotating gas disk disturbed by the jet, gas streaming toward the nucleus, or a ring with another 
smaller counter- or Keplarian-rotating gas disk inside. 

Key words, galaxies: individual (M51, NGC 5194) - galaxies: ISM - galaxies: Seyfert 



1. Introduction 

Active Galactic Nuclei (AGNs) are believed to be powered 
by gas accretion. This gas is supplied from interstellar mat- 
ter in host galaxies, and the gas may form rotationally- 
supported structures around the central supermassive black 
hole. If they are viewed close to edge-on, they may ob- 
scure the central activity from direct view. AGNs can be 
categorized as type 1 if seen face-on, and type 2 if seen 
edge-on; this explanation is known as a unified model (e.g. 
lAntonucci fe Miller! 1 19851 ). Indeed, a few hundred pc res- 
olution molecular gas imaging towar d the central regions 
of the Seyfert 2 gal axies NGC 1 068 dPlanesas et alJll99U; 



ies. A pair of radio jets emanates from the nucleus and 
narro w line regions (NLRs) are a s sociated with the jet 
(e.g., n ' 



Uackson et ail I1993D and M51 (jKohno et all I1996D show 



strong peaks at the nuclei with velocity gradients perpen- 
dicular to radio jets, which suggest the existence of edge- 
on circumnuclear rotating disks. Recent ~ 50 pc resolu- 
tion imaging studies toward NGC 1068 and the Seyfert 
1 galaxy NGC 3227 support this view, s howing more de- 
tailed structures, namely warped disks (jSchinnerer et al.l 
2000a|U). However, observations toward a few low activ- 
ity AGN galaxies with < 100 pc resolution show lopsided, 
weak, or no molecular gas emiss ion toward the nuclei (e.g. 
iGarcia-Burillo et al.l[200l 120051) . 

M51 (NGC 5194) has also been observed in detail 
with molecular lin es in the past, since i t is one of the 
nearest (7.1 Mpc; iTakats fc Vinkdl l2006f) Seyfert galax- 
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, , , Crane fc van der Hulstl Il992t iGrillmair et all 119971: 
iBradlev et all 12004ft . Interferometric images in molecular 
gas show blueshifted emission on the eastern side of the 
Seyfert 2 nucleus, and redshifted gas on the western side 
(jKohno et all 119961 : IScoville et al] |1998f ). This shift is al- 
most perpendicular to the jet axis, and the estimated col- 
umn density is consistent with that estimated from X-ray 
absorption toward the nucleus, suggesting that the molecu- 
lar gas can be a rotating disk and play an important role in 
obscuring the AGN. Interferometric CO (3-2) observations 
suggest a velocity gradie nt along the jet in addit ion to that 
perpendicular to the jet ([Matsushita et aLll2004j ). These re- 
sults imply more complicated features than a simple disk 
structure. We therefore performed sub-arcsecond resolution 
CO(2-l) and CO(1-0) imaging observations of the center of 
M51 to study the distribution and kinematics of the molec- 
ular gas around the AGN in more detail. 

2. Observation and data reduction 

We observed CO(2-l) and CO(1-0) simultaneously toward 
the nuclear region of M51 using the IRAM Plateau de 
Bure Interferometer. The array was in the new A configu- 
ration, whose maximum baseline length extends to 760 m. 
Observations were carried out on February 4th, 2006. The 
system temperatures in DSB at 1 mm were in the range 
200-700 K, except for Antenna 6, for which a new genera- 
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Fig. 1. Channel maps of the CO(2-l) line. The contour lev- 
els are —3, 3, 5, 7, and 9a, where la corresponds to 5.2 mJy 
beam -1 (= 0.96 K). The cross in each map indicates 
the position of the 8.4 GHz radio continuum peak posi- 
tion of R.A. = 13 h 29" 1 52!7101 and Dec. = 47°lT42'/696 
(jHagiwara et al.ll200lb iBradlev et al.ll2004D . The R.A. and 
Dec. offsets are the offsets from the phase tracking center of 
R.A. = 13 h 29 m 52?71 and Dec. = 47° f 1/4276. The synthe- 
sized beam is shown at the bottom-left corner of the first 
channel map. 



tion receiver gave system temperatures of 150-230 K. Those 
in SSB at 3 mm were in the range 140-250 K for Antenna 
6, and 220-550 K for other antennas. Four of the corre- 
lators were configured to cover a 209 MHz (272 km s -1 ) 
bandwidth for the CO(2-l) line, and a 139 MHz (362 km 
s -1 ) bandwidth for the CO(1-0) line. The remaining four 
units of the correlator were configured to cover a 550 MHz 
bandwidth for continuum observations and calibration. The 
strong quasar 0923+392 was used for the bandpass calibra- 
tion, and the quasars 1150+497 and 1418+546 were used 
for the phase and amplitude calibrations. 

The data were calibrated using GILD AS, and were im- 
aged using AIPS. The data were CLEANed with natural 
weighting, and the synthesized beam sizes are 0'.'40 x 0'.'31 
(14 pc x 11 pc) with a position angle (P.A.) of 0° and 
0785 x 0755 (29 pc x 19 pc) with a P.A. of 13° for CO(2-l) 
and CO(T0) images, respectively. Fig. Q]shows the channel 
maps of CO(2-l) emission with a 20.3 km s -1 velocity res- 
olution. The channel maps of CO (1-0) emission show simi- 
lar features to that of CO(2-l) emission with lower spatial 
resolution. Fig. [2] shows integrated intensity and intensity 
weighted mean velocity maps of the CO(2-l) and CO(l- 
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Fig. 2. Integrated intensity (moment 0) and intensity 
weighted velocity (moment 1) maps of the CO(2-l) and 
CO(1-0) lines. The synthesized beams are shown at the 
bottom-left corner of each image. The crosses and the ref- 
erence positions of the R.A. and Dec. offsets are the same 
as in Fig. [TJ (a) The CO(2-l) moment image. The contour 
levels are (1, 3, 5, 7, 9, and 11) x 0.334 Jy beam -1 km s _1 
(= 62.0 K km s" 1 ). (b) The CO (1-0) moment image. The 
contour levels are (1, 3, 5, and 7) x 0.257 Jy beam -1 km 
s- 1 (= 50.6 K km s" 1 ). (c) The CO(2-l) moment 1 image, 
(d) The CO(1-0) moment 1 image. 



0) lines. The noise levels for continuum maps are 1.2 mJy 
beam -1 at 1.3 mm and 0.54 mJy beam -1 at 2.6 mm, respec- 
tively. We did not detect any significant continuum emission 
at either frequency. 



3. Results 

Most of the CO(2-l) emission is detected within ~ 1" 
(34 pc) of the center, and is located mainly on the eastern 
and western sides of the nucleus. There is also weak emis- 
sion located ~ 277 northwest of the nucleus. The overall dis- 
tribut ion and kinematics are consistent with past observa- 
tions (jKohno et al.lll996tlScoville et al.|[T998h . if we degrade 
our image to lower angular resolution; a blueshifted feature 
with the average velocity of ~ 460 km s -1 at the eastern 
side of the nucleus, and a redshifted feature with an average 
velocity of ~ 500 km s -1 at the western side (Figs.[l][2J see 
also Fig. [HJa). W e refer to thes e main structures with the 
same labels as in IScoville et all (11998D (Fig. Hi). 

Our higher resolution images, however, show more com- 
plicated structures and kinematics than the previous low 
angular resolution observations. Molecular gas on the west- 
ern side of the nucleus, SI, is elongated in the north-south 
direction and separated into two main peaks (Sla and b). 
Sla is located 079 (30 pc) northwest of the nucleus, and Sib 
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is l'.'O (34 pc) to the southwest. On the eastern side of the 
nucleus, the molecular gas has an intensity peak 0'.'6 (20 pc) 
to the northeast (labeled S2), which is located closer to the 
nucleus in projected distance than that of Sla/b. 

The feature SI shows a clear velocity gradient along 
the north-south direction, which is shown in Fig. [He) and 
also in the position- velocity (PV) diagram (Fig. [3^,). This 
radient was previousl y suggested by the CO (3-2) data 
Matsush ita et al.ll2004t ) , but the magnitude of the velocity 
gradient is different. The computation of the magnitude of 
the velocity gradient is similar to that used for the CO (3-2) 
data. The fitting result indicates a velocity gradient within 

51 of 2.2 ± 0.3 km s _1 pc -1 , which is larger than that re- 
ported previously, 0.77 ± 0.01 km s _1 pc -1 (the value has 
been modified by the different distance of the galaxy used). 
This difference is partially due to the larger beam size of 
the previous result; the CO(3-2) data set has a beam size 
of 3"9 x T/6 with a P. A. of 146°, and the velocities of S2/C 
and S3 contaminate that of SI. 

The CO(1-0) maps show very similar molecular gas 
distribution and kinematics as those in CO(2-l) maps 
(Fig. [2b, d). Only the western em i ssion was detected in 
previ ous observations (jAalto et al.l 119991 : ISakamoto et al.1 
QUI, but our map clearly shows the emission from both 
side of the nucleus. 

In addition to the previously known features, our CO (2- 
1) image also shows a weak emission near the nucleus with 
a structure elongated in the northeast-southwest direction 
(feature C in Fig. [5^). This structure could be a part of 
S2, since the velocity map (Fig. [2b) and the PV diagram 
(Fig. [3b>) show a smooth velocity gradient, although most 
of the emission in C comes from only one velocity channel 
(419.7 km s _1 map in Fig.[T|). The velocity gradient between 

52 and C is in an opposite sense to that previously seen 
with the lower angular resolution observations mentioned 
above. This structure is not detected in the CO(1-0) line, 
but a hint of a velocity gradient can be seen in Fig. [2H. 

The total CO(2-l) integrated intensity of SI, S2, and C 
is 25.0 1 Jy km s" 1 , and that of SI and S2 in lScoville et al.1 
(|1998D is 33.44 Jy km s" 1 , so th a t our data detected 75% 
of their intensity. IScoville et all (| 19981) detected ~ 50% 
and 20% of the single dish CO(2-l) flux in redshifted and 
blueshifted emission, respectively, so that our data recov- 
ered ~ 25% of the single dish flux. 



4. Discussion 

4.1. Jet-entrained molecular gas 

Our molecular gas data show a clear north-south velocity 
gradient within the feature SI. We suggested from our pre- 
vious study that this velocity gradient ma y be due to molec- 
ular gas entrainment by the radio jet (|Matsushita et afl 
l2004f h Here we revisit this possibility with higher spa- 
tial and velocity resolution data. Fig. 0] shows our CO(2- 
1) image overlaid on the 6 cm radio continuum image 
(| Crane fc van der~ Hulst 1993). The radio continuum image 
shows a compact radio core coincident with the nucleus, 
and the southern jet emanating from there (note that the 
northern jet is located outside our figure). The CO(2-l) 
map clearly shows that SI is aligned almost parallel to the 
jet. In addition, Figs. [2] and [3] show that the velocity gradi- 
ent in SI is also almost parallel to the jet. 




1.0 0.5 0.0 -0.5 -1.0 -1.5 
Dec. Offset [arcsec] 
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R.A. Offset [arcsec] 



Fig. 3. Position-velocity (PV) diagrams of the CO(2-l) 
line. The contour levels are 3, 5, 7, and 9a, where la cor- 
responds to 5.2 mJy beam -1 (= 0.96 K). (a) PV diagram 
along the north-south elongated SI feature (P.A. of the cut 
is 103°). The positions for Sla and Sib are shown with 
labels, (b) PV diagram along R.A. with the cut through 
the Sla, C, and S2 features (P.A. of the cut is 90°). The 
positions for Sla, S2, and C are shown with labels. 



The velocity increases from ~ 480 km s _1 at Sla to 
~ 540 km s _1 south of Sib. This increment is very simi- 
lar to that observed in the NLR clouds along the radio jet; 
Bra dley et al.l (|2004[ ) measured the velocities and velocity 
dispersions of the clouds using the [O III] A5007 line, and 
showed that the velocity of the southern 1" clouds from 
the nucleus are at Vlsr ~ 440 — 590 km s _1 and the veloc- 
ity increases as the clouds move away from the nucleus (see 
Table 2 and Fig. 9 of their paper This velocity range and 
increment are consistent with our data. Furthermore recent 
observations of H2O masers toward the nucleus also show 
a veloci ty gradient alon g the jet with the same sense as our 
results {kagiwara 2007;), in addtion to the good correspon- 
danc e of the velocity range (|Hagiwara et al.ll200lUHagiwaral 
120071 : iMatsushita et all l2004f) . These results suggest that 
the molecular gas in SI (and the NLR clouds and the H 2 
masers) is possibly entrained by the radio jet. These results 
also suggest that some of the material in NLRs is supplied 
from molecular gas close to AGNs. 

Another example of je t-entrained neutral gas is found in 
the radio galaxy 3C293 (|Emont et al.ll2005l ). The velocity 
of H I gas in absorption spectra toward the AGN matches 
that of ionized gas along kpc-scale radio jets. The spatial 
coincidence is not clear, since the spatial resolution of the 
HI data is lower (25'.'3 x 11'.'9) than that of the ionized 
gas data. Our result is therefore the first possible case of 
entrainment of molecular gas by a jet at the scale of ten pc. 

The better resolution of our new CO data allows us to 
revisit the values of the molecular gas mass, momentum, 
and energy of the entrained gas. We derive 6 x 10 5 M Q , 
8 x 10 45 g cm s -1 , and 3 x 10 52 ergs for these quantities. 
These values are about half of the previous values derived 
from the CO(3-2) data, mainly due to the larger beam, but 
the conclusion is similar; the energy of the entrained gas 
could be similar to that of the radio jet (> 6.9 x 10 51 ergs; 



1 We selected the clouds with a ve locity dispersion of le ss than 
100 km s" 1 ; Clouds 3, 4, and 4a in lBradlev et all (|2004 ). If we 
include all the clouds, the velocity is ~ 440 — 690 km s^ 1 with 
a range of velocity dispersion of ~ 25 — 331 km s ; Clouds 2, 
3, 3a, 4, 4a, and 4b. 
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13 29 53.1 52.9 52.7 52.5 52.3 
RIGHT ASCENSION (J2000) 

Fig. 4. The CO(2-l) integrated intensity image (contours) 
overlaid on the VLA 6cm radi o continuum image (greyscale; 
ICrane fc van der Hulstlll993 ). The contour levels, the syn- 
thesized beam, and the cross are the same as in Fig. [TJ 



ICrane fc van der~H ulst 1993), but the momentum is much 
larger than that of the jet (2 x 10 41 g cm s _1 ). One way 
to explain this discrepancy is t hrough a continuous i nput 
of momentum from the jet fsee iMatsushita et alj|20(j4l for 
more detail discussions). 



4.2. Obscuring material around the Seyfert 2 nucleus 

The feature C is located in front of the Seyfert 2 nu- 
cleus, and the CO(2-l) intensity is about 62.0 K km s _1 
(Fig. [2k). Hence the column density can be calculated as 
6.2 x 10 cm~ 2 using a CO-to- H2 conversion factor o f 
1.0 x 10 20 cm~ 2 (K km s" 1 )" 1 ([Matsushita et all I2004D 
and assuming a CO(2-1)/(1-0) ratio of unity. This value is 
far lower than t hat derived from the X-ray absorption of 
5.6 x 10 24 cm -2 (|Fukazawa et al.ll2001[ ). As is mentioned in 
Sectj3l the missing flux of our data is ~ 75%. However, even 
if all of this missing flux contributes to obscuring the nu- 
clear emission, this large column density difference cannot 
be explained. Changing the conversion factor or the ratio 
by an order of magnitude also cannot explain this large 
difference. One way to reconcile this disparity is to assume 
that C is not spatially resovled, in which case the computed 
column density is a lower limit. Alternatively, the obscur- 
ing material preferentially traced by higher-J CO lines or 
denser molecular gas tracers such as HCN may be involved. 
The CO (3-2) intensity in brightness te mperature scale is 
~ 2 t imes stronger than that of CO (1-0) (|Matsushita et all 
12004 ) . and the H CN(l-O) intensity is also relatively stronger 
(HCN/CO ~ 0.4; lKohno et al.lll996[ ) than normal galaxies. 



outer part of a rotating disk as in the AGN unified model. 
However, our images show a more complicated nature, and 
no clear evidence of simple disk characteristics. 

The simplest interpretation is that SI and S2/C are in- 
dependent structures. Since SI is affected by the jet but 
S2/C is not, SI is expected to be located closer to the nu- 
cleus than S2/C, and the projection effect makes the po- 
sition of S2/C closer in our images. Alternatively, S2/C 
may be close to the nucleus, but the entrained gas has 
been already swept away or ionized by the jet. S2/C has 
a velocity gradient, and therefore can be interpreted as a 
streaming gas, presumably inf ailing t oward the nucleus, as 
is observed in the Galactic Center (|Lo fc Claussenl 119831 : 
IHo et al.lH99l . 

SI and S2/C can also be interpreted as a rotating disk 
that is largely disturbed by the jet, and only a part remains. 
According to the velocity gradient along S2/C, the blucsh- 
fited gas is expected at SI, which is the opposite sense to 
the previous suggestion, but the gas shows no signs of it due 
to the jet entrainment. This is possible from the timescale 
point of view; under this interpretation, SI should have 
a blueshifted rotation velocity of ~ 380 km s^ 1 based on 
the velocity gradient in S2/C. SI has a velocity ~ 150 km 
s _1 higher than the expected rotational velocity, and we 
assume that this is the entrained velocity. In this case, it 
takes 2 x 10 5 years to be elongated along the jet by ~ 1" 
or 34 pc. On the other hand, the rotation timescale at this 
radius is about 2 x 10 6 years, an order of magnitude longer 
timescale. The rotating disk can therefore be locally dis- 
turbed by the jet. 

However, the above two explanations have difficulty in 
explaining optical images of the nucleus; the Hubble Space 
Telesco pe images show "X" sh aped dark lanes in front of the 
nucleus (|Grillmair et al.|[T997[ ). suggesting the existence of a 
warped disk or two rings with one tilted far from another. 
An alternative explanation of the dark lanes is that, as 
previously proposed, there is a rotating edge-on ring with 
S2 as blueshifted gas and Sla as redshifted gas. In this case, 
the feature C can be the counterpart of another dark lane, 
which runs northeast-southwest, although C has to be a 
counter-rotating or Keplarian rotating disk to explain the 
opposite sense of the velocity gradient to that of the Sla/S2 
(Sect. [3]). This configuration explains the "X" shape, but 
has a rather complicated configuration, and it is difficult to 
explain why the inner disk C is not disturbed by the jet. 

We imaged the nuclear region of the Seyfert 2 galaxy 
M51 at ~ 10 pc resolution, and we see no clear evidence of 
a circumnuclear rotating molecular gas disk as previously 
suggested. The molecular gas along the radio jet is most 
likely entrained by the jet. The explanations for other gas 
components are speculative, possibly involving a circumnu- 
clear rotating disk or streaming gas. 
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4.3. Molecular gas at ten pc scale from the Seyfert nucleus 

Previous studies suggest that the blue shifted eastern fea- 
ture S2 and the red shifted western feature SI may be the 
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